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ABSTRACT 


The objective of this analytical study was to develop a rapid 
theoretical analysis on approximate half wavelength elements of an 
antenna configured in a co-planar, symmetrical array. A computer pro- 
gram was written employing the method of moments approach to the 
solution of Hallen's integral equation with an approximate two term 
entire domain expansion assumed for the current distributions. With 
the solution of the current distributions on each element, additional 
calculations were made for the input impedance and admittance values, 
field distributions, power gain, and a graphical output of the 


radiation pattern. 
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I. INTRODUCTION 


A. DEVELOPMENT OF THE STUDY 

A goal in the theory of multi-element antenna design is to have a 
tool which easily analyzes the behavior of an array. The equations 
that describe the behavior of a one voltage-driven element with a 
Simple sinusoidal current distribution are readily available in many 
texts [1]. However, upon the introduction of a second or third element 
with different physical dimensions--either voltage driven or parasitic-- 
the behavior without the assumption of a simple sinusoidal current 
distribution becomes rather complex and lends itself to a computer 
oriented solution. 

A large portion of the theory employed in the programming of this 
software package was presented and discussed in a seminar at tne 
University of Mississippi, sponsored by the National Science Foundation 
in 1972 [2]. 

The computer program of this study was written in the XDS Fortran 
IV language. It was specifically to be used in the computer laboratory 
of the electircal engineering department at the Naval Postgraduate 


School, Monterey, California [3]. 


B. ANALYSIS OF THE PROBLEM 


A brief synopsis of the theory of this study is as follows: 


1. The Development of Hallen's Integral Equation 
2. The Solving of Hallen's Equation for the Current Distribution 
by the Use of an Assumed Two Term Approximate Solution Employing the 


Method of Moments 





3. The Development of the Associated Equations for the Relative 
Vertical Plane Pattern, Relative Horizontal Plane Pattern, and the Gain 
of the Array 

4. An Outlining of the Computer Programming Procedures 


C. OUTLINE OF THE RESULTS 

An antenna array of from one to ten elements was modeled in free 
Space. The simple physical parameters of the designed array are the 
input to the computer program. All calculations were made by the 
computer, producing an output of: (a) a listing of the current 
distribution on each element of the array; (b) the input impedance 
and admittance quantities for each of the voltage-driven elements; 
(c) a listing of the relative field distributions about the array; 
(d) a graphics output of the radiation patterns; and (e) the values 


of the gain in magnitude and in decibels. 


D. SCOPE AND LIMITATIONS 

The procedures outlined were comoleted for a symmetrical, co-planar 
array of from one to ten elements in free space witn each element being 
approximately a half wavelength long. No boom effects were included. 
Each of the ten elements may be 11 ог an array of one 
voltage-driven element and nine parasitic (reflectors or directors) 
elements may be configured. The latter includes possible Yagi-Uda 
arrays. The elements may be spaced arbitrarily apart and the center- 
driven voltage may be any complex value. The radii of the elements are 
restricted to the thin wire approximation (a < .01 X). The graphical 
radiation pattern outputs were of the relative vertical and horizontal 


planes. 





Мо optimality procedures were introduced. Further study in this 
area could be extended by introducing a finitely conducting ground 
plane to the model; include elements which are other than co-planar; 
or include elements which are not symmetrical about the axis, but are 


at angles to the axis. 
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11. FORMULATION OF HALLEN'S INTEGRAL EQUATION 

A. THE DIFFERENTIAL EQUATION FOR THE VECTOR POTENTIAL FROM THE 

ЕСИ ЕТЕР ЭЕМ ГУ 

The basis of this study is the calculation of the approximated 
` two-term current in the solution of Hallen's integral equation, as all 
other quantities are determined from this data. In this section, a 
brief discussion of the formulation of Hallen's integral equation is 
presented. The objective is to obtain an expression of the current 
distribution of a cylindrical center-fed voltage driven antenna in 
terms of its length and diameter. 

Consider a voltage driven element as depicted in Figure l. It is 


known that: 


ЕЗ (11-1) 


Е. = E 11-2 
р р ( | 


Е! 
Z 


In this case, the electrical field intensity, E, and the vector 


potential, A, exist only in the Z-component. Therefore, the equation, 





22 
I == (۵ - 143 (11-3) 
becomes, 
ر‎ 
d À 
4 + KR, = јкЕ, (11-4) 
17 m^ 


Equation (II-4) is correct everywhere, where A.(x,y,z) is the vector 
potential and EJ (xy 32) is the electric field intensity. This equation 
is solved on each element of an array of co-planar elements in the YZ 
plane as seen in Figure 2. The radii of the elements in this figure 


are represented as а 
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Evaluating (II-4) on the surface of each of the N elements 


produces: 
d^R (a. y. .z) 2222 s, 
2 nn’ | س2‎ Ё J 7 "DA 
DU ےھ ھت ا ا‎ 
dZ 0 
ш th) (II-5) 


Solving (II-5) with B, and C, being arbitrary constants and Z' being an 


arbitrary point on the 'nth' element gives: 


A, 27272 = C, cos (kZ) + B sin(kz) 


sk f. [E S) sin(Z-S)k dS 
T. a .y,,S) sin q 


0 
П 
0 1+ + + ٦ (11-6) 
simplifying, 
As (a, s, sz) - C cos(kZ) 
укр E. S) sin(Z-S)k ds 
0 J c-p АМЫГ Ё он : 


(= т (11-7) 


Equation (II-7) is achieved with the requirements that the array of 
elements are perfectly symmetric in the YZ plane so that the scalar 
potential is zero at the center of each element. 

Now an expression for the vector potential, A, will be expressed 
in terms of the antenna current and the two expressions will be set 
equal to each other, for the formulation of Hallen's integral equation. 
The only unknowns of the equation will be the current distribution and 


a constant of integration. 





B. THE DIFFERENTIAL EQUATION FOR THE VECTOR POTENTIAL FROM ELEMENT 
CURRENTS 


At an arbitrary point, 'P', in Figure 2, the vector potential can 
be written: 
N 
A(r) = 24 mo (II-8) 
where A. (r) is the vector potential due to the current on the 'ith' 


element and is expressed as: 





- h. -jk|r - ri] 
; Ш 1 1 | E J 1 
z!=-h. jus 
1 1 
where, I 
Е-Е) = |x“ + (y - yi) + (z - z )° |° (11-10) 





In the above equation, x, y, and z are the coordinates of the point 
'P'. The y; and 2: are the coordinates on the 'ith' element. 


Substuting (II-9) and (II-10) into (II-8) produces (II-11). 


n: 


TA кае RI 
2 z!=-h Ui? E E dz. (11-11) 


б-т; 





| 
3j 
N 
ar 
ті 


Equation (II-11) is the vector potential at an arbitrary point in 
space due to all N antenna currents. Evaluating (II-11) along the 


surface of the 'nth' element produces (I1-12). 


N h. 


ЭГ mE j ЕЕЕ 
Alan Yn) = Ga Tq l (zi) e до (102212) 
i=l ziz2-h. Hm i 
1 1 ni 
where: 1 
1217: К? (212 | 
2 С Ё шон - 2) | (11-13) 





C. HALLEN'S INTEGRAL EQUATION 
Setting equation (II-7) and equation (II-12) equal to each other 
along each of the N elements results in Hallen's integral equation 


ЙГ 12). 
5 SDÓKR, LUZ!) 


N 
41 ZJ, т 


ша 


dZ = 


و 


: 7 
Jk : l 
С,со5(К2) + ап Eo (a, yn S) sin(Z-S)k dS 


Mm .ا‎ o N (11-14) 
where: 1 
_ 2 2 е 2 
К 10254 у = E 十 s = yj ) + (Z - Z') | (11-15) 
Since the current in each element is an even function of Z, the 


limits of the integrals can be changed to zero and h;. With this 


change the set of N equations in (II-14) becomes (II-16), where, 


E 
CENE Y c5 Е 
VOU) 1 5111 (К7) (11-17) 
(2) = TT = = (11-18) 
n 30 ` 
ای‎ ЕЕ SRM .)7 7 
Кт) A ыВ 0 ЕЕ ША (11-19) 
ni i R 1; (252°) 7808001 


With this formulation of Hallen's integral equation, it was this 
study's objective to solve the N set of equations (II-16) for the 
unknowns of I, e Ty and the constants of integration C4 .... CN as 


a function of the arbitrary distance along each element Z'. 
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Iii. JHE THEORY OF THE TWO TERM APPROXIMATE SOLUTION TO 
HALLEN'S INTEGRAL EQUATION 





A. TWO TERM APPROXIMATION 

Even though the N set of equations (II-16) may look rather compli- 
cated, there exists only two sets of unknown quantities. Knowing the 
values of the parameters of an array of elements, the currents and 
constants of integration are easily determined. 

The technique of the two-term solution assumes that the current on 


the element may be expressed as (III-1). 
E 
1,(2) = >»), bp Sin Па (h; - 2 P = 2 (111-1) 
al 


Expanding (III-1) produces (III-2). 


1 


1177) = 3 Sin (h; - n t D^ sin cala; (H. - 2 (111-2) 





In equation (III-1), the 'i' represents the number of the element 
and the 'p' is the number of the coefficient in tne sinusoidal expan- 
sion. This expansion sets the number of unknowns at two current 
coefficients and one constant of integration per element. Therefore, 
there are three unknowns per element in the set of equations (II-16). 
Applying (III-2) to (II-16) produces the set of N equations in (III-3). 
For an example, a three element array (N-3) would have nine unknowns 


in the set of three equations. They are, as seen in (III-3), ba, E 


2 pl b3 3 


07, 27 JE 205 


Do % and C 


1 22 


B. METHOD OF MOMENTS 
In order to determine these unkowns, tne method of moments is 


employed [4]. The set of testing functions selected are à set of deita 
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functions, which enables the usage of the point matching procedure. In 
this study, the selected number of matched points per element equal the 
number of unknowns per element, that is, three. The point matching 
process then produces three equations for every one of the N equations 
in (III-3). Each integral in (111-3) must be solved for Z equalling 
one of the designated matched points and Z' existing as the dummy 
variable of integration. The notation used in (III-3) is as follows: 


ni 


h. 
4 (2: 
mp -/ sin |ps/2h. (h. 


I n I I 
تی‎ e (111-4) 
Z'=0 


Here, the 'i' represents the number of the element, 'p' is the 
number of the coefficient in the sinusoidal expansion (in this study 


p=1 or p=2), and 'm' is the number of the matched point on a particular 


element, 'n'. Employing this notation and arranging all the unknowns 

to the left hand side in (III-3), produces (III-5), where there are 
three equations for every one in (III-3). In (III-5) as previously, 

the unknown quantities are the b coefficients and the C, constants of 
integration. Consider again, the example of the three element array 
(N=3) with its nine unknowns. Equation (III-5) becomes nine simultane- 
ous integral equations with nine unknowns. Every 2 becomes a solvable 


integral with the selection of the matched points ( 


C. MATRIX FORMATION 

Rearranging (III-5) into matrix form produces (III-6), which is in 
the form AX=B. The "A" matrix must be filled in with the desired 
number of zeroes, depending on the number of elements, N. The 'A' 
matrix will always be square and be a 3N by 3N size matrix. The 'X' 


and 'B' matrices will always be a 3N by 1 size matrix. Clearly, the 
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'A' and 'B' matrices contain all known quantities and the 'X' matrix 


contains all the unknown quantities. 
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A. ADMITTANCE AND IMPEDANCE 

The solving of the current distribution from the matrix equation 
(111-6) is the essence of this study. From the current distributions, 
the admittance and impedance quantities are found by using the terminal 
current (current at center of element, where the voltage is applied) 


and Ohm's, Law. 


B. RELATIVE VERTICAL AND HORIZONTAL PLANE PATTERNS 
The calculations begin with the solving of the differential 
equation (IV-1) in the far field region [5]. 


j(1) sin e dZ 


dH, 5 ¿SA 


(IV-1) 


The representation of this equation is found in Figure (3). The 
relative vertical plane pattern En in free space is found from (IV-1) 
by equation (IV-2). 


Es = E = ledni, (IV-2) 


The current distribution calculated previously is pictured in 
Figure (4). The approach used in solving (IV-1) was by superposition 


of the two currents. The current Ї may be substituted into (IV-1) as: 


i | ju(t = s/c) 

1, sin [k(L/2 + Z)] e (1У-3) 
and the current 1) as: 

ju(t - s/c) 

I, sin[2k(L/2 + Z)] e (1У-4) 
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FIGURE 4. TWO TERM EXPANSION OF THE CURRENT 
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In Figure (3), the distance 's' to a distant point is found to be: 


S = үр ۔‎ 70506 (IV-5) 


Substituting (IV-3) and (IV-5) into (IV-1) produces an integral 


equation of the form in (IV-6). 


ах 


f e*sinte + bx)dx ээн 


> [a sin(c * bx) - b cos(c * bx)] (IV-6) 
a" + b 


Carrying through the integration and making use of the trigonometric 
identity in (IV-7) and the relationship that h=L/2 from Figure (3), 


sin“ = | - cos^o (IV-7) 


the result is: 


ЖТ» 
Jh cos(khcosa) - cos(kh) 
27 ^ Oar | sino | (1۷-8) 
where * - ME - r/c) 
Relating this to (IV-2) and setting r=1 gives: 
24751 1604 10-02 27004010) 
Бет = 2001 pes — sina | (1V-9) 


Employing (IV-4) and (IV-6) in (IV-1) produces the second components 


of the superposition. 
j 





jl | | 
Hyp Ё = | син. | 05416056). — (1V-10) 
4 - cos 6 
Relating (IV-10) to (IV-2) and setting r=] produces (IV-11). 
Е, : 36011 | Bu | E d (IV-11) 
4 - cos 0 





Now for multiple symmetrical elements, a phase term is introduced due 


to the separation of elements along the Y axis. In taking this into 


2 





account, the time term of I is dropped and with the substitution of 
k=w/c and r = yisins, the form of b becomes that of (IV-12). 


1, = Lp ° (1У-12) 


Finally, substituting (IV-12) into (IV-11) and (IV-9), the far 
field pattern in the Y-Z axis is achieved, using the superposition 
equation (IV-13). 

(IV-13) 


The relative horizontal plane pattern as depicted in Figure (5) is 
found with equation (IV-14). 


к 1 i Jky.cos¢ 
ES = ؤ‎ 1 


à ] 2) 6 (IV-14) 


C. POWER GAIN 
The power gain was calculated for comparison with an isotropic 
source [6]. The gain is equal to the ratio of the power intensity to 


the power density. The gain is then expressed as in (IV-15). 





G = Am (IV-15) 
Win 
where from the Poynting vector, 
2 
e 2 _ 
With, 
NC. р 
Hr = |I| Ro (IV-17) 


The gain with r=] is expressed as in (IV-18). 


G = dE _ (14-18) 


2 
30|1]° R. 
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FIGURE 5. CONFIGURATION OF PARALLEL ELEMENTS 
IN THE XY PLANE 
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V. PROGRAMMING PROCEDURES 


A. ANALYSIS 

With the equations of the unknown available, the writing of the 
software which would produce computer solutions with graphical outputs 
began. The program consists of the main program, one function sub- 
program, and four subroutines. The following is a brief outline of 


the procedures used, in the writing of the program. 


ROUTINE OF STEPS EMPLOYED 

1. The number of arrays to be analyzed, 1 to 99, and the designated 
number of the graphical output device [7] to be used is read in. 

2. The number of elements in the array to be analyzed, 1 to 10, 
and the wavelength in meters is read in. 

3. The applied complex voltage, value of 'Y' along the Y-axis in 
meters, half length of element 'h' in meters, and radius 'a' in meters 
is read in for each element. 


4. Constants are formed including the designation of the position 
ni 
mp 
5. That portion of the 'A' matrix involving the integral S 


of the matched points for the solution of the integral, S 
کرت‎ 
mp is 
solved. The integration routine employs Weddle's Rule [8] which 
approximates the function with a sixth degree polynomial. 


6. The remainder of the 'A' matrix is then completed by 


appropriately filling in the zeroes and the function, 20 


7. The 'B' matrix is completed by solving the function, y 
8. With the 'A' and 'B' matrices formed, the 'X' matrix is solved 


for by using a Gaussian Elimination routine [9] including the 
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decomposing of the 'A' matrix, the solving for the 'X' matrix, and 
then the attempt to improve the 'X' matrix for the most accurate 
solution possible. 

9. With the solution of the 'X' matrix achieved, the values of 
the current are outputed starting at the center (Z=0) and for every 
tenth of the distance along each element until the end of the element 
(Z=h.). 

10. For every element which is voltage-driven, the input impedance 
and admittance values are outputed. 

11, Starting at theta equal to zero degrees, and for every degree 
through 360 degrees, the values of Es] and دم‎ for each element are 
Summed together. 

12. These values are normalized for output and the appropriate XY 
values are calculated for the graphical computer. 

13. Starting at phi equal to zero degrees and for every degree to 
360 degrees, the values of E for each element are summed together. 

14. These values are normalized for output and the appropriate XY 
values are calculated for the graphical computer. 

15. The XY coordinates of the relative patterns to be drawn are 
sent to the graphics computer for display. | 

16. The power gain in magnitude and decibels is computed and 


outputed. 
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VI. CONCLUSIONS 


A. COMPARISON OF RESULTS 
The parameters of two simple antenna configurations were used as 


input data. Results of this study, were then compared for degrees of 


accuracy. 


B. EXAMPLE ONE. TWO SYMMETRICALLY VOLTAGE DRIVEN ELEMENTS 

The first antenna consisted of two identical symmetrically voltage 
driven elements with a wavelength set equal to one meter (300MHz). 
Each element was driven with the complex voltage of 1+j0. The separa- 
tion between the two elements мас .01 meters (y47-0.005, y570.005). 
Both elements were 0.50 meters in length (L-2h). The radius of each 
element was 0.000277 meters. Computer results produced by the program 
of this study for the current distribution are displayed in Figure 6. 
These results show good correlation with those of King's first-order 
theory [10] and a two term sinusoidal approximation as depicted by 
Entler- [11]. 

With this study's calculated current distributions, the input 
admittances, input impedances, relative vertical pattern, relative 
horizontal pattern and power gain in magnitude and decibels were also 


calculated as shown in Table I. 


۶۰۰٠٠٠٠٠٠٦۱۷٦ TWO. SIMPLE SYMMETRIC CENTER FED DUPOEE 
The second example is a simple symmetric center-fed dipole, with a 
wavelength set equal to one meter (300MHz). The complex voltage 


applied was 1430. There is no separation involved in this example 
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REAL PART 


IMAGINARY PART 





=/ бл ЕШ Б 0 205 бай? 7.5 


KING'S FIRST-ORDER THEORY 
cose BUTLER TWO-TERM SINUSOIDAL BASIS 
+ -+--+ RESULTS OF THIS STUDY 


FIGURE 6. CURRENT DISTRIBUTION FOR TWO SYMMETRICALLY 
VOLTAGE DRIVEN ELEMENTS 
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REAL CURRENT IMAGINARY CURRENT DISTANCE ALONG ELEMENT 


(ma. ) (ma. ) FROM CENTER (meters) 
5 28236974 -2 37995399 . 000 
5.16780772 -2.47847437 ‚025 
4.92962568 -2.50660181 ‚050 
4.57696180 -2.45519573 .075 
4.12109759 -2.31881755 ‚100 
3.57492582 -2.09652088 ‚125 
2.95246973 -1.79220618 ‚150 
2.26848157 -1.41511770 .175 
1.53813553 -0.97857729 .200 
0.77681740 -0.50012839 ‚225 
0.00000000 -0.00000000 ‚250 
INPUT ADMITTANCE REAL INPUT ADMITTANCE IMAGINARY 

(milisiemens) (milisiemens) 

5.28236974 -2.37995399 


INPUT IMPEDANCE REAL INPUT IMPEDANCE IMAGINARY 
( ohms ) ( ohms) 


ШЕ 722207827509 70.90930883439 


OMEGA RELATIVE VERT. PATTERN RELATIVE HORZ. PATTERN 
(deg. ) (absolute) (absolute) 
000 1.00000000000 ‚99950656036 
010 „97786160215 ‚99952143819 
020 „01417041942 ‚99956427764 
030 .81632464159 „99962991265 
040 ‚69429231999 ‚99971022788 
050 5580521887] . 99979611278 
060 +3 ШЫ ‚99987663247 
070 .27633644475 „99994227433 
080 211721 00-227 .99998611978 
090 ‚00000000000 . 00000000000 
100 ‚13729063227 ‚99998511978 
110 "2705044475 ‚99994227433 
120 00 le ‚99987663247 
130 23980921997 | ‚99979611278 
140 ‚69429231999 ‚99971022788 
150 „81643564159 „99962991256 
160 ‚91417041942 ‚99956427764 
170 „97786160215 „99952143819 
180 1.00000000000 „99950656036 


The gain of the array is 3.09 = 4.90db. 


TABLE I. COMPUTER OUTPUT FOR TWO SYMMETRICALLY VOLTAGE 
DRIVEN ELEMENTS 
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(y,=0.0). Following the length and radii restrictions, the length of 

the element was 0.25 meters and the radius was 0.0069 meters. 
Computer results produced by the program of this study, are dis- 

played in Figure 7. These results are plotted along with a two term 

Sinusoidal approximation, and a five term sinusoidal approximation 

as displayed by Butler [12] and مگ‎ сс data by Mack [13]. The 


additional data calculated as listed in Example one is shown in Table II. 


р. RECOMMENDATIONS 

More complicated arrays were also analyzed by the computer program 
of this study. These outputs yielded results which were closely 
related to those of other two earn approximations for current 
distributions [14]. 

In Examples one and two, the matched points used were Z=0.0, 
2=0.5h;, and Zzh;. It is recommended in a study by Darko Kajfez [15] 
that for the closest comparable results, using the two term approxima- 
tion of this study, the matched points should be located at 720.0, 

Z-0. 75h. , and =n. 

This study found that the number of increments that the integral 
function is divided into when being approximated by a sixth degree 
polynomial employing Weddle's Rule, will in some cases have an 
unwieldy effect in the solutions of the integrals. It is recommended 
that the integer, 75, be employed in the determination of the number of 
increments. | 

Recommendations for further study are to employ a different complex 
numerical integration routine with the theory of this study. Possibly 
additional programming could be written in order to allow interactive 


graphics for design and testing procedures. 








‘TWO-TERM SINUSOIDAL BASIS 


一 一 一 一 一 一 FIVE-TERM SINUSOIDAL BASIS 
Xxx X X MACK'S MEASURED DATA 
+++ ++ RESULTS OF THIS STUDY 


FIGURE 7. CURRENT DISTRIBUTION FOR SIMPLE SYMMETRIC 
CENTER-FED DIPOLE 
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MAGNITUDE 
(ma ) 


10.069451 76 
10031137742 
1029630974 
.99144960 
251425727 
22725169 
819153077 
.66716441 
.9144935] 
(22742153 
.00000000 


CO — CO On ^4 00 O 0 


-20 
-23 


-30 
= 
-34 


INPUT ADMITTANCE REAL 


(milisiemens) 


9. 40598988 


INPUT IMPEDANCE REAL 


(ohms ) 


9276685977913 


ОМЕСА 
deg.) 


000 
010 
020 
030 
040 
050 
060 
070 
080 
090 
100 
110 
120 
130 
140 
150 
160 
170 
180 


RELATIVE VERT. PATTERN 
(absolute) 


1. 
.97/40371048 
.91249376631 
.81307695801 
‚68974122014 
. 95322994009 
‚41223256076 
22221311. 
. 13498559880 
. 00000000000 
. 13498559880 
2 ۵ 
‚41223256076 
2990322402009 
.68974122014 
.8130/695801 
.91249376631 
.97740371048 

2 


ARGUMENT 
(deg. ) 


.91494045546 
‚88716574162 
-26. 
-28. 


7 
92611960589 


.97158254089 
.68483783763 
.07240413432 
202 
522 
zn 
-36. 


13865883043 
89288284490 
34223309019 
49146393058 


00000000000 


00000000000 


DISTANCE ALONG ELEMENT 


FROM CENTER 


(meters) 


. 000 
1025 
-050 
.075 
. 100 
275 
. 150 
5 
. 200 
“225 
200 


IMPUT ADMITTANCE IMAGINARY 
(milisiemens) 


-3.59460891 


INPUT IMPEDANCE IMAGINARY 
(ohms ) 


35.45193909900 


The gain of the array is 1.64 = 2.15db. 


COMPUTER OUTPUT FOR SIMPLE SYMMETRIC 
CENTER-FED DIPOLE 


TABLE II. 
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RELATIVE HORZ. PATTERN 
(a 


bsolute) 


. 00000000000 
. 00000000000 
.00000090000 
. 00000000000 
. 00000000009 
. 00090000000 
. 00000000000 
. 00000000000 
. 00000000000 
.00000000000 
.00000000000 
.00000090000 
. 00000000000 
. 00000000000 
. 00000000000 
. 00009000000 
. 00000000000 
. 00000000000 
. 00000000000 











АРРЕМПЇХ А 


GUIDE TO OPERATION OF PROGRAM 


Below is a step by step procedure to follow in order to operate the 
computer program, at the Naval Postgraduate School computer laboratory. 
The software requires the use of the digital computer (XDS9300) and 


the graphics computer (AGT). 


A. FORMATION OF THE DECK 

1. The first card is a 'BOOT' card which is found in the card 
Stalls in the laboratory. 

2. The second card is the 'JOB' card which is followed by the 
cards of the program through to the 'DATA' card as shown in the 
enclosed computer program printout. 

3. Following the 'DATA' card is the initial data card which 
contains three digits in the first three columns of the card. The 
first two digits denote the number of arrays to be analyzed (01 to 99). 
The third digit denotes the assigned number of the graphics computer 
to be used (1 or 2). 

4. After the initial data card, the input data of each of the 
arrays is formed. Each array will have one card plus as many cards as 
there are elements in the array. The first card of each array will 
have the number of elements in the array in the first two columns 
(01 to 10). In the next twenty columns, the value of the wavelength 


in meters 1s to be entered. 
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9. Following the 'wagelength' card, there is one card for each of 
the elements. The format for each of these cards is as follows. 


Columns 01 to 15: Real value of complex voltage applied to element. 
Example 1.0 


Columns 16 to 30: Imaginary value of complex voltage applied to 
element. Example -1.0 


Columns 31 to 45: Value of the Y-axis, indicating the position of 
the element in meters. Example 0.0 


Columns 46 to 60: Value of half length of element (h) in meters. 
Example 0.25 


Columns 61 to 75: Value of radius of element in meters. 
Example 0.001 


6. AS an example, in order to run example one and two of section 


VI of this study, the cards after the 'DATA' card would appear as 


follows: 

Бара |: 021 

ог 

Card 1: 022 

Wand 2: O 120 

and 39: OTO -- 75 
Card 4: EO 020. 072025 7 
Вака 5: I IE 

Card 6: 1202020000020, 25 9069 


B. PREPARATION OF XDS9300 DIGITAL COMPUTER 
1, Press IDLE switch at the XDS9300 control panel. 
2. Hold the RESET switch depressed and press the POWER switch in 


order to energize the computer. 


3. Turn switch at teletype to the ‘on’ position. 


4, Press the READY button on the line printer. 
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C. PREPARATION OF THE AGT GRAPHICS COMPUTER 

1. Press ON at the operator's control panel (OCP). 

2. Wait 5 seconds, then press RESET on the OCP. 

3. Press the ON/START switch at the disk drive of the AGT to be 
sed, 

4. Press HALT, RESET, and BTSP, in that order on the OCP. 

9. Verify that the 'location counter' (behind the door of the 
left-most rack) indicates zero or all lights out. If not, press 
RESET and BTSP again. 

6. Load a paper tape marked DISK BOOT LOADER (located behind the 
door of the left-most rack) into the tape reader of the AGT's teletype. 

7. Verify that the knob on the front of the AGT's teletype is in 
the 'line' position. 

8. Put the tape reader switch to 'start.' 

9. Now verify that the 'location counter' indicates 110 (octal). 
If it is not, return to step 1 of this section and begin the 
procedure again. 

10. When step 9 is satisfied, verify that the switches in the 
center of the OCP have the configuration 24000 77776 (octal). 

11. Press HALT, RESET, RUN, and PULSE 1 on the OCP. 

12. If the AGT's teletype does not respond now with a request 
for the operator to input MO/DA/YR, press HALT, RESET, RUN, BUTTON 29, 
and PULSE 1. If teletype does not respond now, return to step | of 
this section. 

13. Type RESET("GATED",101)! on the AGT's teletype. 

14. Type GATED: on the AGT's teletype. 


15. Turn vertical and horizontal gain knobs fully clockwise. 
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D. READING IN PROGRAM DECK 


1. 
2, 


Place cards in card hopper. 
Push POWER ON button, at the card reader. 
Push START button, at the card reader. 


Push IDLE, RESET, RUN, and CARDS on the XDS9300 control panel. 
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